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O . ABSTRACT 

■ In a systematic search over 11 cluster fields from Cluster Lensing And Supernova 

survey with Hubble (CLASH) we identify ten passively evolving massive galaxies at 
redshift z ^2. We derive the stellar properties of these galaxies using HST WFC3/ACS 
multiband data, together with Spitzer IRAC observations. We also deduce the optical 

O \ rest-frame effective radius of these high redshift objects. The derived stellar masses and 

measured effective radii have been corrected by the lensing magnification factors, which 
c/^ I are estimated by simply adopting the spherical NFW model for the foreground cluster 

^ lens. The observed near-IR images, obtained by HST WFC3 camera with high spatial 

resolution and lensed by the foreground clusters, enable us to study the structures of 
such systems. Nine out of ten galaxies have on average three times smaller effective 
^ . radius than local ETGs of similar stellar masses, in agreement with previous works 

■ at redshift 1.5 < z < 2.5. Combined with literature data for z ~ 2, we find that the 
\ mass- normalized effective radius scales with redshift as Ve/M^-^^ oc (1 + z)~^-^^. We 

confirm that their size distribution shows a large scatter: from normal size to 5 
times smaller compared to local ETGs with similar stellar masses. The 1-a scatter 
^logve of the size distribution is 0.22 and 0.34 at 2; ~ 1.6 and z ~ 2.1,respectively.The 
observed large size scatter has to be carefully taken into account in galaxy evolution 
CN ■ model predictions. 
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1 INTRODUCTION to be much smaller, by a factor of - 3 - 4 (Daddi et al. 

, ^ , _ ^ , , . . 1 2005; Trujillo et al. 2006; Toft et al. 2007; van der Wei et 

Local Early-Type Galaxies (ETGs) are spheroidal systems ^^ ^008; van Dokkum et al. 2008, 2009; Damjanov et al. 2009; 

which predominantly consist of old stars with mass-weighted g^^^^^^ 2009; Cassata et al. 2010, 2011; Ryan et al. 

ages 01 > 8 — 9 Gyr see Henzini 2006 tor a review). Inis or^o at ^ 1 or^o o ^ 1 or^o rv- 

. . ~ .11 r 1 . T^rr^^ 2012; Newman et al. 2012; bzomoru et al. 2012; Zirm et 

provides strong evidence that most 01 the stars m ETGs , or^oN rrom /TT^T^r-io • 1 ^1 ^ t-.t^^ 

^ „ 1 1 1 .p r 1 11 . 1 al. 2012). Deep HST/WFC3 images reveal that PEGs at 

were formed at redshitt z > l.b and evolved passively. Re- ^ ^ iii^ir-^^ii 1 

, . T 1 ^r^r^oN r f 1 . z ^ 2 are not surrounded by the tamt extended envelopes 

cent observations e.g. Kriek et al. 2008) tound red, massive r ^ • w o ^ 1 or^ox ^ xi 1 j 

^ ^ o ^ T PI T^rr^^ 1 • 01 material (e.g. bzomoru et al. 2012) . On the other hand, 

ETGs already at z > 2. Many or these ETGs are desig- ^1 1 i r ^ • 1 • ^1 1 1 -r^ • i- 

^ . ^ ^ . - . /T^-^^ X 111 tl^6 lack ot compact massive galaxies at low redshiit implies 

nated as passively evolving galaxies (PEGs), generally by a - -n ^ • 1 ^- r o ^ /rn --n ^ 1 

^ . . r j/ i -(^ f significant size evolution from 2 ~ 2 to z ~ (Trujillo et al. 



2009). 



color-color selection criterion and/or a low specific star for- 
mation rate (sSFR). Massive PEGs contribute a significant 
fraction to massive galaxies at 2; > 2, with a number density Two popular models for the size evolution have been 

of 1 X 10""^ Mpc~^ (Kriek et al. 2008). proposed: (1) dissipationless (dry) minor merger (Naab et 

Sizes of local ETGs have been found to tightly correlate al. 2009, Hopkins et al. 2010); (2) "puff-up" due to the gas 
with stellar masses (Shen et al. 2003). While compared to mass loss by AGN (Fan et al. 2008, 2010) or supernova feed- 
local ETGs of similar stellar masses, the sizes of massive back (Damjanov et al. 2009). However, recent studies found 
(M^ > 10^°-^ M0), PEGs at redshift z ^ 2 have been found that observed size evolution of the compact massive galax- 
ies at redshift z > 1.5 can not be explained well by the 
present models. Nipoti et al (2012) constructed a ACDM- 
^ E-mail:llfan@ustc. edu.cn based analytic framework, supported by suites of N-body 
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merger simulations, to predict the evolution of ETGs un- 
dergoing dry mergers. Using a compilation of observations 
of ETGs at 2; = 1 — 2.5, they concluded that mergers alone 
are not consistent with the observed rate of structural evo- 
lution at 2; > 1.5. For the puff- up model, Ragone-Figueroa 
& Granato (2011) found the mass loss by AGN feedback can 
produce a significant size increase based on the results of nu- 
merical simulations. However, they augued most of the puff- 
ing up occurs when the stellar population are much younger 
than the estimated ages of compact high-z PEGs, so that 
the puffing up due to the gas mass loss may not dominate 
the size evolution observed so far. 

Size evolution with redshift has been described by re oc 
(1 -h z)^ in the literature (e.g. Newman et al. 2010, Dam- 
janov et al. 2011, Cimatti et al. 2012). However, the value 
of P is diverse for different works, varying from ^ —0.8 to 
—1.6. The large diversity of /3 may be due to the differ- 
ent sample selection criteria (morphological, color-color or 
sSFR selections) and/or partly due to the sparsity of data 
points at high redshift. For instance, there are only about 
ten objects at 2; > 2 in Damjanov et al. 2012. Using UV J 
and sSFR selection criteria, Szomoru et al. (2012) selected 
21 PEGs at ^ - 2 in the GOODS-South field as part of the 
CANDELS survey. They found those 21 PEGs at ^ - 2 not 
only had smaller effective radii than local ETGs of similar 
stellar masses, but also had a significant spread in their size 
distribution (criogre ^ 0.24). If the larger size scatter is real 
at higher redshift, it will be another observational constraint 
for dissecting the size evolution models. 

In this letter, we will use multiband data from the 
CLASH survey to search massive PEGs at 2; 2. 
HST/WFC3 F160W-band images will be used to derive the 
sizes of PEGs. And complementary with the recent archive 
data, we will study size distribution at z ^ 2 and size evo- 
lution with redshift. Throughout this letter, we assume a 
concordance ACDM cosmology with Q^n = 0.3, Ha = 0.7, 
Ho = 70 km s"^Mpc"\ All magnitudes are in the AB sys- 
tems. 



2 DATA AND SAMPLE SELECTION 

The Cluster Lensing And Supernova survey with Hubble 
(CLASH) is an HST Multi-Cycle Treasury program that is 
acquiring images in 16 broad bands with the Wide Field 
Camera 3 (WFC3) UVIS and IR cameras, and the Ad- 
vanced Camera for Surveys (ACS) for 25 clusters (Post- 
man et al. 2012). The images are reduced and mosaiced 
(0."065 pixel" ^) using the MosaicDrizzle pipeline (Koeke- 
moer et al.2011). Photometric redshifts are estimated for 
all the galaxies using the full 16-band photometry via both 
the BPZ (Benitez 2000) and LePhare (Arnouts et al. 1999) 
codes. In this letter, we use the released photometric cata- 
log of 11 clusters Q- We supplement the HST observations 
with Spitzer IRAC data. We use a software package with 
object template-fitting method (TFIT; Laidler et al. 2007) 
to measure IRAC photometry. TFIT uses the spatial posi- 
tions and morphologies of objects in a high-resolution image 
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Figure 1. Broad band SEDs (black points) and best-fit models 
(solid lines) of z ~ 2 PEGs. On the top-left of each panel, the 
number marks the galaxy ID in our sample (see Table 1). 



to construct templates, which are then fitted to a lower reso- 
lution image. The fit can be simultaneously done for several 
close objects in a crowded field. We use the F160W band 
as the TFIT high-resolution template to measure the IRAC 
photometry. 

Stellar masses and star formation rates (SFRs) are esti- 
mated from SED fitting to the full photometric data set (see 
Figure 1), assuming a Kroupa IMF (Kroupa 2001) and the 
BC03 stellar population models (Bruzual & Chariot 2003). 
We use a Bayesian-like code GIG ALE (Noll et al. 2009) to 
carry out the SED fitting. As described in Noll et al. 2009 
, an exponentially decreasing star formation history and a 
modified Calzetti attenuation law have been adopted (see 
Calzetti et al. 2000 for more details). The SED-derived stel- 
lar masses and SFRs are corrected by the foreground lens 
magnification factor. We simply adopt the spherical NEW 
(Navarro et al. 1997) model for the foreground cluster lens. 
The lensing properties of NEW halo can be calculated an- 
alytically for given position (Bartelmann 1996). The mass 
and concentration of the lens clusters are obtained from the 
literature (Umetsu et al. 2012, Zitrin et al. 2011, Gavazzi et 
al. 2003 and Kling et al. 2005). The detailed properties of 
the lens clusters are not taken into account, such as mass 
ellipticity, substructures, etc. The estimated magnifications 
may be slightly underestimated, especially for the clusters 
not long after merging. 

We select passively evolving galaxies using a selection 
criterion based on sSER, with sSFR< 0.3/ti/, where tn is 
the Hubble time. As additional selection criteria, we also 
demand that the PEGs will be at 1.5 < z < 2.5 and stellar 
masses > 10^^'^ Mq. The final sample consists of ten 
galaxies. Our selection criterion is very close to those using 
a fixed sSER (e.g. sSFR < IQ-^^yr'^ ). And we also check 
if our selection criterion is consistent with the color-color 
selection. We find eight out of ten PEGs in our sample fulfill 
the J — 3.6/im and z — J color-color selection (Papovich et al. 
2012). J — 3.6/im and z — J colors correspond approximately 
to rest-frame V — J vs. U — V at z ^ 2, which Williams et 
al. (2009) showed can effectively separate passively evolving 
galaxies from star-forming galaxies. 



^ http:/ /archive. stsci.edu/prepds/clash/ 
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Table 1. The Sample of Passively Evolving Lensed Galaxies at 1.5 < z < 2.5 



Cluster^") 


id(^) 


Ra^'^) 


Dec^'^) 


0(d) 






log M^^^ 


log sSFR(^) 


n(0 




b/a^'^) 






[deg] 


[deg] 


[arcsec] 






[Mo] 


[yr-^] 




[kpc] 




MACSJ1206-08 


1 


181.53618 


-8.81180 


64.7 


1.83 


1.998 


10.85 


-11.15 


6.3 


1.4 


0.39 


MACSJ0647+70 


2 


101.99016 


70.26104 


57.4 


1.13 


1.556 


11.08 


-10.88 


3.9 


5.1 


0.75 


MACSJ0647+70 


3 


101.99335 


70.23792 


55.6 


1.15 


1.941 


10.99 


-11.23 


3.1 


0.9 


0.50 


MACSJ0647+70 


4 


101.99383 


70.23749 


57.3 


1.15 


1.942 


11.20 


-11.17 


6.0 


1.5 


0.53 


MACSJ0744+39 


5 


116.21017 


39.43837 


74.0 


1.14 


2.276 


10.88 


-10.91 


1.8 


0.7 


0.62 


MS 2137-2353 


6 


325.07814 


-23.67539 


71.2 


1.38 


1.534 


10.74 


-11.16 


2.2 


1.0 


0.36 


RXJ1347-1145 


7 


206.90139 


-11.74505 


84.9 


2.03 


1.646 


10.77 


-10.76 


2.0 


2.2 


0.64 


RXJ1347-1145 


8 


206.86622 


-11.76696 


64.1 


3.30 


1.811 


10.56 


-11.27 


4.5 


2.0 


0.32 


RXJ1347-1145 


9 


206.88078 


-11.77017 


58.3 


3.84 


1.662 


10.54 


-10.96 


3.0 


1.3 


0.47 


MACSJ1149+22 


10 


177.40348 


22.41853 


74.4 


1.09 


2.007 


11.04 


-11.13 


4.0(0 


0.9 


0.87 



(a): The foreground cluster; (b): ID of PEGs; (c): Right ascension and Declination of PEGs; (d): Projected angular distance of PEG from 
the cluster center; (e): The lensing magnification factor which is estimated by simply adopting the spherical NEW profile for the foreground 
cluster lens. The mass and concentration of the lens clusters are obtained from literature (Umetsu et al. 2012, Zitrin et al. 2011, Gavazzi 
et al. 2003 and Kling et al. 2005); (f): Photometric redshift of PEGs using BPZ ; (g): Stellar masses , which have been homogenized to 
BC03 models and Kroupa IMF ; (h): Specific star formation rate ; (i): Sersic index ; (j) : Effective radius, which has been corrected by the 
lensing magnification factor ; (k) : Axis ratio ; (1) : Sersic index is fixed to 4.0 . 
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Figure 2. HST F160W band original (left) and residual (right) 
images of passively evolving galaxies in our sample. The size of 
each panel is 6 x 6 arcsec^. On the top-right of each panel, the 
number marks the galaxy ID in our sample. 



3 SIZES OF PEGS AT Z - 2 

We use the GALFIT package (Peng et al. 2002) to fit Sersic 
profiles to the PEGs in the HST F160W band. The fits are 
performed in a 6 x 6 arcsec^ region (See Figure 2). Empiri- 
cal PSFs instead of model PSFs are used in the fits. We ob- 
tain the empirical PSF for each cluster using PSFEx (Bertin 
2011). PEG-1 is close to another galaxy. Two independent 
Sersic models are simultaneously fitted to PEG-1 and its 
companion, respectively. We use the same procedure to fit 
PEG-3 and PEG-4 (both of them are PEGs). The Sersic 
index n is fixed to 4 for PEG- 10 in order to make the fit 
converge. 

The observed images of high redshift PEGs in our sam- 
ple have been lensed by the foreground clusters. We need to 
take the lens magnification into account in order to estimate 
the intrinsic sizes of PEGs. For simplicity, we assume that 
the lensed image has the same light profile as the source. 
Following the definition of the magnification factor /i, we 
have 



where Ts is the intrinsic angular size of the source image, is 
the measured angular size of the lensed image, both in unit of 
arcsec. /i is the lensing magnification factor estimated using 
the method introduced in the Section 2 (see Bartelmann 
& Schneider 2001 for more detail on lensing). The measured 
angular size from the best-fit Sersic model is corrected using 
the formula in Equation 1 and then converted to physical 
size using photometric redshift. 



4 DISCUSSION 

van der Wei et al. (2011) argued the majority of PEGs at 
2; 2 in their sample were disk-dominated. We investigate 
this argument with our sample. Half of our PEGs have best- 
fitting axis ratios b/a < 0.5 and six (PEG-1, PEG-3, PEG-4, 
PEG-6, PEG-8 and PEG-9) have b/a < 0.6 (See Table 1). 
The axis ratio distribution indicates that a significant frac- 
tion of PEGs at z ^ 2 may be disk-dominated. However, 
we find that Sersic index n are mostly greater than 2, indi- 
cating elliptical-like morphology. Only two galaxies (PEG-5 
and PEG-7) have n < 2. We re-fit our PEGs with two Sersic 
models. We fix the index n of one Sersic model to 4 to rep- 
resent the bulge component and n = 1 to represent the 
disk component. Two galaxies (PEG-7 and PEG-9) show a 
significant disk component with > 50% light in disk. The 
results by axis ratios and two-component model show that 
20 — 60% of PEGs at z ^ 2 are disk-dominated in our 
sample. 

In Figure 3, we plot the distribution of effective radius 
versus stellar mass for our PEGs sample at 1.5 < z < 2.5. 
In our sample, both stellar mass and effective radius are 
corrected by the lensing magnification factor. We also com- 
pare with the local relation (Shen et al. 2003) and other 
literature results with similar range of redshift and stellar 
masses (Cimatti et al. 2008, Cassata et al. 2010, Szomoru 
et al. 2012, Ryan et al. 2012 and Zirm et al. 2012). All stel- 
lar masses are homogenized to BC03 and Kroupa IMF. We 
find that our results are consistent with literature results at 
z ^ 2. Nine out of ten PEGs in our sample are below the 
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Figure 3. The distribution of effective radius versus stellar mass 
for passively evolving galaxies at 1.5 < z < 2.5. Our PEGs sam- 
ple is represented by black filled circles. Both stellar masses and 
effective radius are corrected by the lensing magnification fac- 
tor. The zblcr errors of stellar masses and effective radius are 
given by GIG ALE and GALFIT, respectively. Other literature 
data are over-plotted with plus signs (Gassata et al. 2010), aster- 
isks (Gimatti et al. 2008), diamonds (Zirm et al. 2012), triangles 
(Szomoru et al. 2012) and squares (Ryan et al. 2012), respec- 
tively. All of the stellar masses are homogenized to BG03 and 
Kroupa IMF. The solid line shows the local effective radius and 
stellar mass relation of early- type galaxies (Shen et al. 2003). The 
dashed lines indicate the ±lcr scatter of Shen et al. (2003) rela- 
tion. The arrows on top-left show the moving directions due to 
the uncertainty of the lensing magnification factors fi. 



local relation, with re 1 kpc. These galaxies are on aver- 
age '-^3 — 4 times smaller than the local counterparts with 
similar stellar masses. The compact sizes of high redshift 
PEGs can be formed naturally according to the dissipative 
collapse of baryons (see Fan et al. 2010). 

In Figure 3, we find that the sizes of PEGs at ^ 2 
span a wide range and need a significant evolution to catch 
up the local relation. Furthermore, we explore the size evo- 
lution with redshift combining our sample together with the 
literature data. We divide the combined data into two red- 
shift bins, z ^ 1.6 and z ^ 2.1. The former redshift bin in- 
cludes 25 PEGs and the later bin has 31. The average stellar 
masses of both bins are 8 x 10^° M©. We adopt the local 
size-mass relation re oc M^, where a equals 0.56 by Shen et 
al. (2003). We plot the average mass- normalized sizes of two 
redshift bins in Figure 4 (filled squares). We also plot the 
local relation with shaded area in Figure 4. We use the rela- 
tion re oc (1 + z)^ to describe the redshift evolution. For the 
combined data, we have /3 ^ — 1.13±0.13 (See the dashed 
line in Figure 4) . The /3 value is in agreement with the result 

of Gimatti et al. (2012). They derived /3 1.06 ± 0.14 for 

the mass- normalized radius with logM^/M© > 10.5. How- 
ever, other results with steeper slope had been found from 
different samples. Using optical and near- infrared photom- 
etry in the UKIRT Ultra Deep Survey and GOODS-South 
fields of the GANDELS survey, Newman et al. (2012) found 
the size evolution of quiescent galaxies with redshift could be 
described as \og{re/M^-^^) = 0.38-0.26(^-1). This relation 
can be approximately rewritten as re/M^'^^ oc (1 + 2;)^, with 
P ^ —1.33. Damjanov et al. (2011) found an even steeper 
value, 13 ^ —1.62 ±0.34. The sizes need to increase by a fac- 



Figure 4. Evolution of the effective radius with redshift. The 
shaded area refiects the distribution of local SDSS galaxies. The 
filled squares show the average sizes of PEGs at z ~ 1.6 and z ~ 
2.1. The error bars show the la scatter. The dashed line represents 
the size evolution with redshift described by re oc (1 + z)"-'^ -'^'^. 
The other symbols are the same as those in Figure 3. 



tor of - 3 for /3 - -1.0 from z = 2 to 0. While for /S 1.6, 

the sizes would increase by a factor of ^ 6 during the same 
comic epoch. The diversity of observed size evolution with 
redshift is at least partly due to the sparsity of data points 
at 2; ~ 2. A large number sample of PEGs ai z ^ 2 will be 
helpful to deicide the valid size evolution with redshift and 
constrain the different galaxy evolution models. 

A relevant feature of the data for massive PEGs at high 
redshift is the quite large scatter of the size, as is apparent 
from Figures 3 and 4. We plot the la size scatter (aiogre) of 
two redshift bins in Figure 4. We find that the size scatter is 
large at 1.5 < z < 2.5 in the combined data. The values of 
c^iogre are 0.22 and 0.34 at 2; ^ 1.6 and z ^ 2.1, respectively. 
Newman et al. (2012) showed a very similar result laiogre — 
0.24 and 0.26 at 1.0 < ^ < 1.5 and 1.5 < ^ < 2.0, respec- 
tively.In Figure 4, the error bars show la size scatter of two 
redshift bins and the shaded area reflects the distribution of 
local SDSS galaxies. Gompared to the small size scatter of 
SDSS ETGs (cTiogre = 0.16 at 2; ^ 0.06), the size scatter also 
needs a significant evolution with redshift (about a factor of 
1.5 - 2 from z - 2). 

Larger size scatter at higher redshift need to be ac- 
counted for by galaxy evolution models. If the observed 
size evolution with redshift is dominated by the dry mi- 
nor merger, we will expect that more compact PEGs have 
higher dry merger rates relative to less compact objects at 
given redshift range. Another possibility is size evolution at 
2; ^ 2 could be dominated by AGN feedback. One assump- 
tion is that gas will be removed by the strong AGN feedback 
in a short timescale. In this case, galaxy structure after AGN 
feedback will relax to the new equilibrium. Another situation 
for AGN feedback mechanism is that the gas removal will not 
complete in only one single shot. Martizzi et al. (2012) used 
cosmological simulations to show the effect of AGN feed- 
back on the gas distribution in the central regions of cluster 
galaxies. They found gas heated and expelled by AGN feed- 
back can return after cooling; repeated cycles could modify 
the stars and dark matter mass profile. During the cycles, 
the observed sizes will oscillate from compact (all gas is lo- 
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cated within the central regions of galaxies) to normal (gas 
is expelled from the galaxy center). This produces the large 
size scatter observed at high redshift, and the size scatter 
will be related to AGN activities which decrease from z 2 
to 2; 0. For AGN feedback mechanism, larger gas fraction 
at higher redshift will also be needed in order to explain the 
observed size scatter evolution. Most recently, Olsen et al. 
(2012) found the ubiquitous presence of AGN in massive, 
quiescent z ^ 2 galaxies. This result provided observational 
indications for the important role AGN played in explaining 
the size (and size scatter) evolution. 

It should be noticed that the lensing magnification cor- 
rection that we adopted is simple, which may introduce un- 
certainties as well. However, the uncertainty on the lens- 
ing magnification may not significantly alter our main re- 
sults above. The corrected stellar mass is proportional to 
while the corrected effective radius is proportional to 
l/y/71. As jjL decreases/increases, the location on the size- 
mass plane will move along the directions shown by arrows 
in Figure 3. The arrows (ve oc M^'^) nearly parallel to the 
local size-mass relation (re oc M°"^^). So the uncertainty of 
II will have a small effect on the mass-normalized effective 
radius. For instance, the mass-normalized effective radius 
with ^ 10^^ M0 will only change by a factor of < 5% 
when /i increases by a factor of two. 



5 SUMMARY 

In this letter, we use the CLASH survey data to select a 
massive, lensed PEGs sample at 1.5 < z < 2.5, based on 
the SED-fitting sSFR and accurate photometric redshift. We 
deduce the optical rest-frame effective radius of these high 
redshift objects using GALFIT. Combined with literature 
data, we find the size evolution of passively evolving galaxies 
scales with redshift as re oc (l+z)-^-^^^°-^^ We confirm that 
their sizes have a large spread. The 1-a scatter aiogre of the 
size distribution is 0.22 and 0.34 at 2; ^ 1.6 and z ^ 2.1, 
which is about 1.4 and 2.1 times larger than the local value 
(c^iogre 0.16), respectively. The observed large size scatter 
has to be carefully taken into account in galaxy evolution 
model predictions. 
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